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Abstract 

We study the jj^ expansion in the recently proposed model of the lattice gauge theory 
induced by heavy scalar field in adjoint representation. In the first approximation the 
fluctuations of the density of eigenvalues of the scalar field are Gaussian, so that the 
scalar glueball spectrum is defined from the corresponding linear wave equation. 







1 Introduction 



The problem of solving large QCD seemed hopeless until recently some new approach was 
taken0, ^. Instead of trying to solve it starting from the asymptotic freedom, we suggested 
to start much earlier, at such scales, that there are no gluons yet, but rather some gluon 
constituents. 

The idea is that one could then use the freedom of choosing any constituents and any 
form of their interaction, as long as the effective theory at large scales would be still QCD. 
Clearly the gauge symmetry should be there from the start, but there need not be the Yang- 
Mills (or Wilson, in the lattice formulation) term in the bare Action. This term should be 
induced in an asymptotically free region, as result of integrating out heavy constituents. 

This idea is quite old. Perhaps Sakharov was the first one to suggest the induced Gravity. 
After that the idea was unsuccessfully implemented by the author (see the review paper [Q) in 
the fermionic string model of QCD. The heavy constituents were introduced as 2-dimensional 
Majorana fermions at the world sheet of string. Formally, the Schwinger-Dyson equations 
of the large QCD were satisfied, but one could not make any sense out of this model. 

As we see it now, after ten more years of study of the string theory, the real problem 
comes from the string. Apparently, there is no such thing as a noncritical string in any 
number D > 1 of external dimensions. The interaction at the world sheet must be essentially 
nonlocal, and/or it should involve some extrinsic geometry. The idea of string realization 
of the Lorentz group was a beautiful one, and it still has chance to work in Gravity, but in 
QCD I personally do not believe in this idea any more. Apparently, the coordinates could 
only be arguments of fields, rather then fields in QCD. 

Coming back to the induced QCD, the idea of Kazakov and myself was to take heavy 
scalar field in adjoint representation of the gauge group as a constituent field. Naively, if 
one just integrates this field out, one induces the correct Yang-Mills action, with positive 
coefficient in front, scaling correctly at large A^. By the way, the simpler choice of the 
fundamental representation would yield too large induced gauge coupling, so it could not 
induce an asymptotically free QCD at large N. The large (~ A^) number of fundamental 
scalars would do it, but this model is hard to solve. 

The correct computation of induced gauge coupling is not so trivial, as there are feedback 
effects from induced hard gluons, and those from effective scalar quartic interaction. If 
one starts from asymptotically free region, and goes to smaller scales, there would be two 
conflicting effects. First, there would be antiscreening of gauge coupling from the gluons, 
which would tend to decrease the gauge coupling further. Second, there would be screening 
from the scalars, which is numerically smaller, then antiscreening, so that alone it cannot stop 
the logarithmic decrease of effective gauge coupling. The desired infinite gauge coupling at 
lattice scales could be generated by scalar quartic interaction, which is known to grow without 
limits at small scales. At least we know, that asymptotic freedom is not consistent with 
scalar quartic interactions. At larger scales, the scalars' decouple, being heavy, so that true 
asymptotic freedom is left intact. So, the scalars effectively cut off the asymptotic freedom 
at the scale of their mass, and move the effective coupling towards another "confinement" 
region at lattice scales. 
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What happens in this model at small scales is beyond the reach of the perturbative 
field theory. However, if we take the particular nonperturbative theory, we have a chance 
to adjust the parameters so, as to visit the perturbative region on the way from small to 
physical scales. In fact, it is hard to imagine how could we miss the asymptotic freedom, if 
we would get nontrivial theory in four dimensions. The common belief is that there is only 
one gauge theory in four dimensions: the asymptotically free, quark confining QCD. 

The unique property of the scalar field in adjoint representation is that corresponding 
lattice gauge model without the Wilson term is soluble at large A^, for arbitrary scalar 
potential. The one link integrals over the gauge fields reduce at large to certain 
functional determinant, which was analyzed in previous works0, ^. In particular, in the 
second work[^, we found some coupled integral equations, exactly soluble by means of the 
Riemann-Hilbert method. 

The gauge field integration produces effective action, which depends only on the density 
of eigenvalues of the scalar field. One could therefore, integrate out the unitary matrices, 
which diagonalize the scalar field, or simply take the gauge where this field is diagonal. The 
corresponding Faddeev-Popov Jacobian is nothing but the Vandermonde determinant, which 
produces in effective Action the term, quadratic in the density of eigenvalues. 

At A^ = oo this density does not fluctuate: it is given by the solution of the classical 
equations, minimizing this effective Action. By some historical reasons it was difficult for 
everybody, authors included, to accept this simple and natural scenario. We were misled 
by large A^ matrix models of 2D Quantum Gravity, where the Itzykson-Zuber determinant 
was used in a different context. There, it reduced to the Slater determinant, thus revealing 
the hidden fermionic degrees of freedom. In higher dimensions there is no such thing as 
fermionization of the Bosonic problem, neither could there be a spatially inhomogeneous 
vacuum. So, there seems to be no alternative to the space independent vacuum density of 
the eigenvalues of the scalar fleld, which we call the Master Field. 

Still, the technical problem of reexamining the old C = 1 matrix model with the new 
representation of the Itzykson-Zuber determinant was pending, so the author (unpublished) 
checked within the strong coupling expansion, that the old fermionic solution can also be 
correctly reproduced by the spatially independent Master Field. Recently, this problem was 
completely solved^], so that there is no longer any contradiction between the fermions and 
the Master Field. 

The classical solutions for the density come in two phases: the strong coupling phase, 
with density continuous at the origin, and the weak coupling phase, where there is a flnite 
gap in the spectrum. The critical phenomena near the transition point are much richer 
then in the usual one- matrix models. The critical indices, computed in|Q are transcendental 
numbers at D > 2, in particular, the density grows as p(A) ~ A" and physical mass scale 
777,2 r-u where cosvra = — 3^^ , a > 1. 

The scaling laws in the four-dimensional theory are amazing, but they do not contradict 
the asymptotic freedom, as there is extra renormalization from the lattice scales to the 
region of the asymptotic freedom. In fact, the induced gauge coupling Pq = is expected 
to come out as some negative number times logarithm of the scalar mass in the lattice units, 
according to the KG analysis at the small scale region. Comparing this with the usual RG 
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analysis in the large scale region, we, indeed, get the scaling laws for the physical mass in 
terms of the bare scalar mass. 

There is, however, one unsolved problem in the induced QCD, that of the spurious local 
Abelian symmetry. It is well known, that the lattice theory in adjoint representation has 
extra local Zn symmetry. This issue was analyzed a long time ago|0, ||, and it was noted, 
that unbroken symmetry locally confines quarks so that they cannot move separately even 
at the lattice scales. It was suggested in the recent work^, that this symmetry breaks 
spontaneously at the transition point, as it does in the mean field approximation]^, |[. In 
any case, this is the problem for the higher approximations, as the center of the gauge 
group is negligible at = cxd. 

The present paper is devoted to the next correction in the scalar sector of the strong 
coupling phase of the induced QCD. 

In Section 2 we develop the general method of expansion, and find exact integral 
equations for the Gaussian fluctuations of the density. 

In Section 3 we go to the local limit, using the same method, as before, and we find 
another Riemann-Hilbert problem, soluble in the scaling limit. 

In the Appendix we present the computation of the second variation of the Itzykson-Zuber 
determinant. 



2 ^ Expansion 

At infinite there are no fluctuations of the density of eigenvalues of the scalar field, and 
therefore, no interesting physics. These fluctuations show up in the first order of the 
expansion. There are two methods for the systematic expansion: the Schwinger-Dyson 
equations and the saddle point method. 

For our present purposes the saddle point method appears to work better, so we elaborate 
it here, using the modification suggested by G.Parisi (unpublished). The idea is to change 
the integration variables from the eigenvalues to the corresponding density 

M^) = ^E^(A-A,) (1) 

j 

The change of variables goes as follows 

J d^Xoc J VpJ[p] (2) 

where 

J[p] = J Veexp (^tN J dXe{X)p{X)^ J d^Xexp (3) 

We neglect everywhere the constant normalization factors, as we are interested only in av- 
erages. 
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It is convenient to shift e(A) by ^lnp(A) , and then single out the translational mode 
e = const which yields the density normalization condition 



6[1-J dXpiX) j (4) 
The translational mode can be eliminated by the background gauge condition 

6 ( [ dXp{X)e{X)) (5) 



with the trivial constant Jacobian, which does not require ghosts. After simple transforma- 
tions we find 

J[p] = 5(1- j dXp{X)^ exp (^-N J dXp{X) Inp(A)^ J[p] (6) 

where 

J[p] = J Ve6 (I rfAp(A)e(A)) exp (-^ J dXp{X)e'{X)^ (7) 

exp N(^tJ dXe{X)p{X) + dXp{X)e\X) + In (^J dXp{X) (exp(-ze(A))))) 

Now we could expand the last exponent in e, keeping in mind, that density fluctuations 
p — p ^ jf and e ~ A^^2. The resulting expansion of J is straightforward, this is one 
dimensional perturbation theory with propagator 

In the leading order we could set J = 1 + 0[jj:], which leaves us with extra term 

N f dXp{X) In p{X) (9) 



in effective Action. This term is local and linear in density, so it effectively shifts the scalar 
field potential, without altering the quadratic form of the second variations. 

This quadratic form can be computed exactly, using the technique of the previous work[^]. 
Let us outline this computation. First let us rederive the saddle point equation with the new 
technique. The normalization condition would be satisfied identically for 

p(A)=,(A)-^;^ (10) 

where il^{X) vanishes at the endpoints of the support of p{X). The first variation with respect 
to tp{X) is equivalent to differentiation ^ and setting Aj A, as one can readily check. 
The effective Action for the x - dependent density pr^. reads 

Seffip] = (11) 

-N'Y. J dXp,{X) J dX'p,{X')\n\X-X'\ + N'Yl J dXp,{X)U{X) - ln/[p.,P,] 



<xy> 



T.{nJ d\px{\) lnp(A)+lnJ 
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where the first term comes from the Vandermonde determinant, the second one from the 
scalar field potential, the third one from the Itzykson-Zuber determinant, and the remaining 
terms - from the above Jacobian. These last terms would contribute only in the higher orders 
of the expansion. 

The background density of eigenvalues Q is determined at = oo from the saddle point 
equation 

^Seff[p] _ 1 d 6Seff[p] 



5^(A) NdX 6p{X) 
This equation is the same, as before 



(12) 



m - -^^'^''.y'^^ (13) 

where 

V'{X) = [ dX'-^ (14) 



A- A' 

and the function F{X) represents the logarithmic derivative of the Itzykson-Zuber deter- 
minant. This function for arbitrary density satisfies the set of Schwinger-Dyson equations, 
which were derived and studied in the previous work. Substituting the classical formula for 
F{X) in these equations, we arrived]^ at the Master Field Equation (MFE) 

ft(A) = ^C/'(A) + ^^»r(A) (15) 

»nA)^/fa.c.a„^ (16) 

Our objective now is to compute the second variation of the effective Action. The poten- 
tial term does not contribute, being linear in density, the Vandermonde determinant yields 
a local term 

-pj dXm J dX'^iX')^^-^. (17) 
The Itzykson-Zuber determinant yields two terms, the local one being 

-D J dX^{X) J rfAV(A')a(A, A') (18) 

with 

and the nonlocal one 

- ^ / dXij^iX) I dX' J2 V'.+^(A')^(A, A') (20) 



^ In the rest of the paper we denote the background density as p rather then p. 
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with 

Collecting all the finite terms, and going to the local limit, we find the standard form of 
effective action of continuum field theory 

S^ff = J d^x£{x) (22) 

with effective Lagrangean 

£{x) = JdxJ dX' ^ri{X, A') d^iP{x, A) A') (23) 

^ ^ +D (a(A, A') + r7(A, X'))] ^(x. A) ^(x. A') 



+ 



{X'-xr 

,P'(A) 



J dXilj{x,X)'- 



P(A) 

The linear term shifts the vacuum average of the density by 



MA) = -^V'^(A) (24) 



where ipo satisfies the equation 



The particle spectrum is described by the wave equation for the fiuctuations if} — ipo — 

Y{X)e'^^ 



O^JdX' Y{X') ( + ^^(A, A') +{D- ip^)r/(A, A')) 



(26) 



In the Appendix we study the perturbed Schwinger-Dyson equation for the second vari- 
ations and we find the following linear integral equations for the kernels a, rj 

jd\'K(X„, X'HX', A") = (27) 

/<iA';C(A„.A').(A',A") ^ / .A'K(Ao.AO (-^^ . (28, 
where ^^^^ 

^'^ ^ (Ao-i?(>/))^ + 7rV(A0 ^^^^ 

and 

''<^«- - A^» (/ ^(A^ 1^^) 



Now we eliminate the unknown functions a, r] by integrating the wave equation over A 
with the weight /C(Ao, A) which yields 



^A^Ao,A)^ 

^^-2^n jm^^ = y y ^) ( (A^ - GGVA) iv j 

(32) 

These equations should be solved together with the MFE. 

All above equations were exact in a sense that we did not go to the local limit, except for 
replacing the lattice derivatives by local ones f\, neither did we assume anything about the 
support of eigenvalues. The same equations hold for the weak coupling phase, where there 
is a gap at the origin. In this paper we study the strong coupling phase where there is no 
such gap. 



3 Local Wave Equation 

Let us consider the local limit, when the density at the origin vanishes. The analysis of 
the Master Field Equation shows, that the scaling solution, independent of the ultraviolet 
cutoff, requires, that p(A) and r(A) = R{X) — X vanish as some power of the cutoff, so that 

p{u) ~ r{u) < u (33) 

In this case the MFE can be expanded in p, r as follows 

/oo p^iy) 
dv — (34) 
-oo (A — uy 

r{X) = ^u'{X) + ^^V'{X) 

where renormalized potential u'{i>) starts from the linear mass term m^u. The higher order 
terms are irrelevant in the scaling region, but they might be important for the full MFE, 
to provide required cancellations of the regular terms in the scaling region. These regular 
terms represent substruction terms in dispersion relations. There is no need for these terms 
in the full MFE at the lattice, but in the local theory there are ultraviolet divergencies. 

This convenient method of elimination of regular terms is to introduce two analytic 
functions 

= 20D) - '^^) 

^One could restore the lattice theory by replacing by 2^^sin^ ^P^. 
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/oo 
(i/i (36) 
-oo IX — Z 

with the symmetry property 

V{-z) = -Viz); Qi~z) = -Q{z) (37) 
and note that at '^z +0 by construction 

<^Q = {'^Vf. (38) 
On the other hand, in virtue of the above equation for density 

3? Q = ^-^'^ (^') (39) 

This is the nonhnear Riemann-Hilbert problem. 

The key identity, used in derivation of local MFE is the following one 

rfA'/C(A, X')AiX') AiX) - (r(A)A(A))' + p [ "^^^^^l^f^ + Oip') (40) 



(A - A') 

where the first two terms come from the small region A' — A ~ p(A), and the last one comes 
from the region A' ~ A. The integral in the small region can be reduced to the residue at 
the complex pole of the kernel. To be more precise, the O(p^) correction to this formula was 
also used in derivation of the local MFE, where the linear terms cancel, but we do not need 
this correction below. 

The expansion of the function G{X, 0) up to linear terms in p reads 

GM)^ ,^ } .J i+^i ''T'r' . ) (41) 



(X-4'-r(4>)) V 1 (4,-i')(\-i>)^ 



Substituting this into the wave equation (|3^), we find variety of terms. To reduce them, it 
is convenient to introduce the analytic function 

such that 

QH<P + ^0) = M^)p f^P^ (43) 



The function Y{(p) can be reconstructed from dispersion relation 



FW = /d. ^-^^^+^'°> ln|^-.| (45) 
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Let us now turn to the wave equation. After some algebra, we reduce it to the following 
boundary problem^ 

2^p(0)3?^'(0) + (-F^ + 2 + 2D^ V"{(t))) 53^(0) + 2{D - l)p(0) J^lM^M^ = q (46) 

Let us now solve these problems in the scaling limit. The local MFE allows the scaling 
solution, in proper units 

7rp(A) = A°; a = 1 H — arccos — — (47) 

TT 31/ — 2 

Various branches of the arccosine correspond to various fixed points of the theory. This 
scaling form holds for A ^ rrC where 7 = and m is the physical mass scale. 
In this case the functions ^'(A), r(A) are given by 

^V\X)^ I^A- A"tan-7ra (48) 

r(A) ^ ^ ^™ (49) 

The relation between the mass m\ and the scalar potential is as follows 

ml = U"{0) = 2D + u"{0) = 2D + 2 - + (50) 
so that the critical point is at 

^ + D-1 
mlc = (51) 

In this paper we shall find the solution at ^ mj. Then we are left with the scaling 
terms, so that the scaling Ansatz for JF could be used 

J^{z) = I'i-izY; s = {0,1} (52) 

The discrete parameter s corresponds to the parity 

J^i--z) = i-lYHz) (53) 

and the ratio of real and imaginary parts follows from symmetry and analyticity. Substituting 
this Ansatz into the equation, we find the equation for the index e 

= . tau (^) H- ^^.) + . cot (54) 

There are infinitely many solutions Cg^n which grow asymptotically linearly with the principal 
quantum number n. 



The real and imaginary parts at real are understood, as usual limits from the upper half plane. 



9 



In general, there are also powerlike corrections to the density, so that the solution is not 
so elementary. The complete computation of the spectrum requires an exact solution of the 
nonlinear Riemann-Hilbert problem for the density. This is the subject of the next paper. 

One important thing is clear right away: this is not a free field theory at D > 1, neither 
we see any pathology at D > 1, like those of the Liuoville theory. This is confining theory, 
but not quite the string. In general, there are infinitely many masses, and they grow with 
quantum numbers, as one would expect for QCD. 

4 Discussion 

The puzzle still remains unsolved. Is this QCD? 

At least this is the theory of scalar particles, confined by strong interaction with gluon 
field. Would gluon field decouple, like it docs in one dimension, we would get trivial free 
particle spectrum. Instead, we are going to get infinitely many particles, with growing 
masses, stable at infinite N. We expect the decays to show up in the next ^ approximations. 

In QCD, we would like to get the mesons, from confined qq pair, in addition to pure 
glueballs. However, the quarks are locally confined in this model, unless the Zjy symmetry 
would break, spontaneously or otherwise. This seems to be the most urgent thing to do. 

Also, we would like to see the perturbative QCD. This can be checked already in the 
scalar sector. If there is admixture of glueballs to the scalar branch of the spectrum, then the 
correlation function of fiuctuations of the scalar eigenvalue density should have logarithmic 
singularity at oo, times some power of P^. This singularity in perturbative QCD comes 

from the two gluon exchange, via the operator. Here it must come out as something like 

T—, ^- InP' (55) 

P^ + const n ^ ^ 

due to the divergencies of the spectral sums. 

Everything seems to be set up for the comparison of induced QCD with perturbative 
QCD. Let us work hard and keep the fingers crossed. 
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A Second Variation of Itzykson- Zuber Integral 

In the previous paper we derived the following coupled set of linear integral equations for 
functions Gx{p),F{v) Q 



^'^"^ - ^'^""^ = -1 + (A - F{v))GM (56) 
jj, — V 

d^p,{^)GM = I dv^^ (57) 
J X — V 



where F{p) is the first variation of the logarithm of the Itzykson-Zuber integral (see the 
section 2.). The densities of eigenvalues at two endpoints of the link were denoted as pi and 
p2. In order to find the second variations, we should vary this equation with respect to pi 
and p2. 

Let us compute nonlocal variation rj first, as it is technically simpler. In this Appendix 
we denote the eigenvalues as 0, p and A. We are going to vary F, G with respect to p2- These 
variations 

d 6F{<P) d 6Gx{<P) 

dil)6p2{ilj) dip dp2{ip) 

satisfy the differentiated version of above equations ([56|), (|57|) 



Viij, u)G,iu) = (A - Fiu))g,ii/j, u) - [ dpp(^)^^i^liil_M^ (59) 

dpp{p)g\{ip,p) = (60) 
As in0, we have to introduce the auxiliary analytic function 

= ,61) 



In this Appendix we treat A as subscript rather as argument of G. 
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in the upper half plane of z and we find from above equations the following linear relations 
for the corresponding boundary values 

Ax^^{z^i^) {\- R{z) -iTrp{z))-Ax,biz-iO) {X - R{z) + in p{z)) = -27np{z)Gx{z)r]{i), z) 

(62) 

At infinity this function decreases as j with coefficient in front known from (|5D|) 

MA^) - I (63) 

The solution of the boundary problem is straightforward. We use the previous solution 
of the homogeneous problem 

Gxip) = Y^^f ^a(^ + ^0) (64) 



Tx{iy + iO) _ A - R{u) + mp{u) 
Tx{u-tO) ~ X - R{u) - i7ip{u) 



(65) 



Tx[z) = exp / — -arctan- — - (66) 

' ' 7r(z/ — z) A — R[v) ) 



and we find 



Ax,M = Tx{z)Bx,^{z) (67) 

where the new function Bx,xi){z) which decreases at infinity precisely as Ax,xi){z) is to be 
reconstructed from imaginary part 

c-R ^ ^ -TTp{z)Gx{z)'q{:ip,z) 

^ + «0 = — — 68 

Tx[z + zO) (A - R[z) - ircp{z)) 

This expression can be simplified in virtue of above equations as follows 

is Bx,i,{z + zO = — , — 69 

[X — R[z))^ + 7r'^p^{z) 

so that the resulting Cauchy integral reads 

= /d.^*) ^ (70) 

J [X — R[z)Y + 'K'^p'^[z) z — V 

Finally, comparing the asymptotic behavior of this integral with known asymptotics 
7^(oo) = 1 and required condition ( pBD for Ax,^{z) we find linear integral equation for 

P(0)??(^,0) _ 1 ^^^-^ 



" (A-i?((/)))2 + 7rV(0) (A-^)2 
Let us now turn to the second cj) derivative, and differentiate the above equations for the 
first derivative. We shall use the same notation gx{<p',(f)) for derivatives of the G function 
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The differentiated equations read 

^ (73) 
dlJ'PifJ')9xi(p',fJ') = (74) 

This equation differs from the previous one only by the right side. Repeating the same 
steps as before we find 

z-^ - ^^^"^ y '^'^ 7^ (A-i?(0))2 + ,rV(0) ^ ^ 

Finally, setting 2 oo in above formula and using the equation (|7^) we find equation 
for a{(j)' , (f)) 

p(0)(a(0^0) + ^4.-^M5fM) 



(A-/2(0))2 + 7rV('/') 



(76) 
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